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Abstract 



Matsumoto and Yoshimura ||hep-ph/9910393| have argued that there 



are loop corrections to the number density of heavy particles (in thermal 
equilibrium with a gas of light particles) that are not Boltzmann sup- 
pressed by a factor of e~^^l'^ at temperatures T well below the mass M of 
the heavy particle. We argue, however, that their definition of the num- 
ber density does not correspond to a quantity that could be measured 
in a realistic experiment. We consider a model where the heavy parti- 
cles carry a conserved U(l) charge, and the light particles do not. The 
fluctuations of the net charge in a given volume then provide a measure 
of the total number of heavy particles in that same volume. We show 
that these charge fluctuations are Boltzmann suppressed (to all orders in 
perturbation theory). Therefore, we argue, the number density of heavy 
particles is also Boltzmann suppressed. 
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In a series of papers, Matsumoto and Yoshimura (hereafter MY) have challenged the 
conventional wisdom concerning the number density of a gas of heavy particles in thermal 
equilibrium with a gas of light (or massless) particles. They consider a model of a heavy 
spin-zero boson (represented by a real scalar field (f) interacting with a massless spin-zero 
boson (represented by a real scalar filed x)- The scalar potential is 

Vi^, X) = \M'v' + + ^A,x' + \X^\\ (1) 

where the couplings are all real and positive, and we take A<^ ^ A^, A ^ A^, and A^ < 1. 
This hierarchy among the couplings allows the x particles to function as an efficient heat 
bath for the (p particles. 

At the level of free field theory, the equilibrium number density of (p particles is given by 

where g = 1 counts the number of species of (f particles, 

/(^) = ^ (3) 

is the Bose distribution function, /? = 1/T is the inverse temperature, and E = (p^-|-M^)^/^ 
is the single-particle energy. For T <^ M, we have 

n^o = {2^)-"\MTfl\-P''- (4) 

the factor of e~^^ means that n^o is Boltzmann suppressed. However, MY argue that there 
are loop corrections to n^o that are not Boltzmann suppressed; specifically, they find 

= n^o + cX^T^M^ + ... (5) 

for T -C M, where c = 1/69120, and the ellipses stand for all higher-order corrections. 

They key issue that we wish to address (raised also in [Q) is the underlying definition 
of n^. For T -C M, MY define via = p^/M, where is the energy density of the ip 
particles. This energy density is in turn defined (for all temperatures) via 

where H is the total hamiltonian, and 

= + liViff + |mV^ + counterterms (7) 

is the free-field part of the ip hamiltonian, plus counterterms (some of which involve the x 
field) that are necessary to remove infinities in this composite operator. 

Eq. (|^) is a highly plausible definition of p^. However, it does not correspond in any 
obvious way to how the number density of ip particles would be determined experimentally. 
Standard methods all involve a search for individual, on-shell ip particles. Real-world ex- 
amples of this include present-day dark matter searches, and measurements of the cosmic 
microwave background radiation. 
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What is needed theoretically is a measurable attribute that is carried by the particles 
only. To create one, we modify the model slightly by making the Lp field complex, and 
requiring its interactions to conserve the corresponding U(l) charge. We leave the x fi^ld 
real, and the x particles neutral. The modified scalar potential is 

V{v, x) = M'if^if + iA^(y.ty,)2 + ^A^x' + I^V^VX'- (8) 

We can now study the net charge contained in a large but finite volume V. Of course, in 
thermal equilibrium, the average net charge (Q) vanishes, but it has nonzero fiuctuations 
(Q^). If we weakly gauge the U(l) symmetry with a small (and therefore dynamically 
irrelevant) gauge coupling e <^ A, we can in principle measure these charge fiuctuations 
without tracking individual ip particles. 

For T -C M, it is easy to compute (Q^)o, where the subscript indicates that we are (for 
now) neglecting interactions. The number iV+ of positively charged particles in a volume V is 
then controlled by a Poisson distribution; this implies (A^+)o — (^+)o — (^+)o- The number 

of negatively charged particles is controlled by an independent Poisson distribution, with 
{N^)o — {N-)o = {N-)o. Overall charge neutrality implies (A^-)o = (^+)o = ^n^poV, where 
n^o is now given by Eq. (^ with g = 2. The net charge is Q = — N_, and so 

{Q')o = {{N^-N^Y)o 

= (iV2)o+(iV2)o-2(iV+)o(iV_)o 

= n^oV. (9) 

We see that the charge fiuctuations give us a measurement of the total number of particles 
in a given volume. At higher temperatures (but still ignoring interactions), quantum effects 
modify this result to 

The ratio {Q'^)o/n^oV depends weakly on temperature; it rises slowly from one at T <^ M 
to 7r^/6C(3) = 1.37 at T ^ M. (The result is similar for fermions, with a ratio of one 
at low temperatures and 7r^/9C(3) = 0.91 at high temperatures.) Furthermore, it seems 
highly unlikely that weak interactions could significantly modify Eq. (^. If we were to have 
either (Q^) -C n^V or (Q^) ^ n^V, we would be forced to conclude that the movements 
of positive and negative particles are highly correlated (in order to suppress or enhance the 
charge fiuctuations in any particular volume). This is inconsistent with the usual notion of 
a gas of particles that move freely and independently between occasional scatterings, and 
would appear to require strong interactions. 

We therefore propose to define the number density of particles, for T <^ M, via 

= miV. (11) 

This definition has the advantage (not shared by the definition used by MY) of being di- 
rectly connected to the experimentally measurable quantity (Q^). Adopting Eq. (|TI|) as 
our definition of n^, the question becomes whether or not the loop corrections to (Q^) are 
Boltzmann suppressed. 
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To compute these loop corrections, we introduce a chemical potential /i and the partition 
function 



Z = Tie-^^"-^'^\ (12) 

We then use 



i2 



(13) 

^l=o 



At the one-loop level (that is, for free ip particles, and ignoring the /^-independent contribu- 
tions of the X particles), we have the textbook formula 

\nZo = vJ ^ln[l + /(E-/i)] + (/i^-/.). (14) 

Using Eq. ([ID in Eq. (|l|) yields Eq. (0). 

As in , loop corrections to In Z may be computed via finite temperature perturbation 
theory (see, e.g., 0). At the two- and three- loop level, the contributing diagrams are shown 
in Fig. (1). Vertices joining four ip lines are neglected since we have assumed that ^ A^. 
Furthermore we neglect diagrams with only x lines, since these are independent of fi. The 
(p propagator is 

where n is an integer specifying the discrete energy. (With M = /i = 0, this is also the x 
propagator.) The direction of the flow of charge is the same as the flow of energy. A closed 
(p loop with a single vertex contributes a factor of 

T A{n,p) = —[l + f{E + f,) + fiE-f,)], (16) 

n=-oo ^ ^ 

where E = (p^ + M^Y^^. The "1" in square brackets yields a divergent, temperature 
independent term that must be removed by renormalization. The remaining two terms are 
Boltzmann suppressed, since at low temperatures 

f{E ± /i) ~ e-^*'^e^^^e"^P'/2*^. (17) 

This implies that the two-loop diagram in Fig. (1), and the first two three-loop diagrams, 
are all Boltzmann suppressed. The third three-loop diagram has a factor of J2n ^i^^ P)^) 
and turns out to be suppressed as well. This leaves only the last ("basketball") diagram. 

According to the standard rules of finite temperature perturbation theory the con- 
tribution of the basketball diagram is 

— InZ,, = CA T Ey ^ ^ (2.) 6 (p, + p. + ks + k,) 

n s 

x/35„,+„2+„3+„4,oAiA2A3A4. (18) 
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FIG. 1. Two- and three-loop contributions to In Z/f3V that include propagators for both the 
heavy f particle (solid lines) and the massless x particle (dashed lines). 
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Here pi and p2 are the momenta on the if hnes, and ks and k4 are the momenta on the x 
hnes. Each momentum and energy is taken to flow from the left vertex to the right vertex. 
Thus we have, in the individual propagators, masses Mi = M2 = M and M3 = M4 = 0, and 
chemical potentials /xi = —^2 = and 11^ = ^4^ = 0. Note especially that /i2 = — /ii; this is 
because, with our convention on energy flow, the charge flow must be opposite to the energy 
flow on one of the lines, and this changes the sign of the chemical potential. Finally, 



C 



1 
2! 



■ 2 



1 
4 



(19) 



is a combinatoric factor. The 1/2! comes from the expansion of exp(— /3ifint), the (—1/2)^ 
comes from two vertex factors arising from the last term in the scalar potential, Eq. (^, 
and the 2 comes from the number of ways to match up the x lines. 

Eq. (|18D can be evaluated by the standard procedure p of writing the Kroneker delta 
as an integral. 



/35ni+...+n4,0 



g27rj(ni+...+n4)t//3 



(20) 



and then performing each sum via contour integration. 



E 



(27rn//3-i/i)2 + E2 2tt 

= JL 

2E 



— (h dz 



1 



(21) 



Here the z contour originally encloses the real axis, and then is deformed into two small 
circles surrounding the poles aX. z = ± E); this deformation is allowed for < t < /5. 
Putting all of this together we now have 



In Zbb = CX^ I dpi dk-s dh i2n)^S%pi + pa + kg + k4) / dt D1D2D3D4, (22) 



where dp^ = d^pi/ (27!-)^ {2Ei), and 



^(l3-t)(E,+tj.i) 



7(^^ + /i.) + e*(^'-^')/(^.-/iO 
The integral over t is tedious but straightforward to perform; the result is 

H^i + /i+)('5i + /2-)(53 + ms, + u) 



(23) 



j dt D1D2D3D4 =iiiiii 



5i=0 53=0(54=0 



SiEi + £1-^2 + ^3-^3 + £^4-E'4 



+ 



(5l + /i+)(l - 5i + /2+)(53 + + h) 



+ (/i -/i). 



(24) 



where Ei = (—1) % and /14 
repeated use of the relation 



SiEi — £1-^2 + ^sE^ + e4^E4 
f{Ei + fi), /2_ = f{E2-fi), /s = /(^s), etc. We have made 



e^^fiE) = l + f{E) 



(25) 
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in obtaining Eq. (0). 

We are interested only in those terms (if any) in Eq. ( ^4|) that are not Boltzmann sup- 
pressed. As we see from Eq. (p!7D, any term containing a factor of either fi± or f2± is 
Boltzmann suppressed. Thus, we can drop all such terms. At this point, we see immediately 
that the remaining terms are all independent of fi, since the only /i dependence in Eq. ( p^ is 
that which is contained implicitly in the factors of fi± and /2±. Therefore, all /^-dependent 
contributions (at the two- and three-loop level) to In Z are Boltzmann suppressed. Eqs. (|T3p 
and ([TI|) then imply that (Q^) and are also Boltzmann suppressed. This is our main 
result. 

For completeness, and for comparison with the results of MY , we continue the com- 
putation of the unsuppressed terms in InZbb- We set fi± = f2± = 0. All numerators in 
Eq. (^) are then either 1, /a, f^, or 73/4. The first of these yields a divergent, temperature 
independent term that is canceled by a renormalization of the vacuum energy. The second 
and third yield divergent, temperature dependent terms that are canceled by renormaliza- 
tion of the (f and x thermal self-energies (see the corresponding discussion for QED in 0). 
Dropping these terms, we have 

■^0 e3=±l e4=±l J^l + -^2+ ^Z^-i + ^4-t^4 

Substituting this into Eq. (|2^) , making a low temperature expansion, and dropping a final 
divergent contribution to the vacuum energy ultimately yields 

J_ InZbb = , + 0(Ti°). (27) 

13V 648000M4 ^ ^ > ^ ' 

The corresponding correction to the total energy density p is then given by 

This is consistent with the results of MY who also find that the leading correction to p 
is of order T^/M^. 

That there are corrections of this form to p is not surprising. Without interactions, 
we have p = (7r^/30)T^ from the x particles, plus the Boltzmann-suppressed contribution 
of the ip particles. At temperatures T <^ M, we should be able to integrate the heavy ip 
field out of the functional integral, and be left with an effective lagrangian for the massless 
X field alone. This lagrangian will contain nonrenormalizable interaction terms that are 
suppressed by powers of A and inverse powers of M. These will give rise to corrections like 
Eq. (pSl). From this point of view, it is clear that all such corrections should be thought of 
as modifications of the energy of the x particles, and not as unsuppressed contributions to 
the energy of the particles. 

We now argue that (Q^) (and hence n^p) is Boltzmann suppressed to all orders in per- 
turbation theory. Consider an exact evaluation of the partition function 

Z = J2e-^^^--f'^-\ (29) 
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where the sum is over a basis of energy and charge eigenstates. States that yield a /i- 
dependent contribution to Z must have Q 7^ 0. States consisting of a single stable ^ 
particle are well defined exact energy eigenstates, with E = (p^ + M^)^/^ and Q = ±1; their 
contribution to Z is obviously Boltzmann suppressed. Furthermore, any other state with 
Q = ±1 must have energy E > M; otherwise, the ip particle would not be stable (since it 
could decay into this lighter charged state). At tree level, energy eigenstates with two or 
more if particles all have E > 2M; interactions can modify this to ii^ > 2M[1 + 0(A^)]. The 
O(A^) term could in principle be negative (if there is a bound state of two or more like-charge 
particles), but cannot make E <C 2M in a weakly coupled theory. We therefore conclude 
that all //-dependent contributions to Z are Boltzmann suppressed. We have already seen 
this explicitly at the level of two and three loops. 

With this in mind, we can return to the original model of MY, Eq. (|1]). At the level 
of Feynman diagrams, the models are essentially identical, the only difference being the 
combinatoric factor associated with each diagram. We therefore expect that should be 
Boltzmann suppressed in this model as well. We are still in need, however, of a general defini- 
tion of that does not rely on the trick of Eq. (pJ]), but which corresponds to experimental 
measurements. 

We conclude, in accord with [^], that the proper definition of the number density of 
heavy particles (in equilibrium with a thermal bath of light particles) is a delicate matter. We 
have argued that it is essential that be defined in a manner that renders it measurable in a 
realistic experiment. In this paper, we have considered a model in which the heavy particles 
carry a conserved U(l) charge, while the light particles are neutral. In this situation, the 
local charge fluctuations provide an experimentally accessible measure of the number density 
of heavy particles. We have shown explicitly that, at low temperatures, this number density 
is Boltzmann suppressed up through three-loop order in perturbation theory, and we have 
argued that this must in fact be true to all orders. 
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